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I.  INTRODUCTION 


Most  transient,  one- dimensional,  interior  ballistic  calculations 
have  assumed  that,  once  ignited,  the  propellant  regression  rate  is  given 
by  a quasi-steady  power  law  dependence  on  pressure.  Increasing  the  low 
pressure  rate  by  adding  an  arbitrary  constant  (of  the  same  order  as  the 
quasi-steady  rate)  has  occasionally  improved  the  agreement  between  pre- 
dicted and  observed  pressure  waves  in  guns*^’  . Such  low  pressure  en- 
hancement is  typical  of  the  transient  burning  rates  which  have  been  pre- 
dicted to  exceed  the  quasi-steady  rate  by  several  multiples^  for  a KTSS 
type  model  and  by  about  one  multiple^  for  the  Zeldovich  model.  This 
paper  will  report  use  of  the  Zeldovich  transient  burning  model  in  the 
NOVA  flamespreading  computer  model  applied  to  a 155mm  howitzer  and  to  a 
five  inch  gun  to  predict  an  effect  on  the  ignition  and  combustion  of  a 
granular  bed. 

The  Zeldovich  approach  assumes  a quasi-steady  gas  phase  in  which  the 
transient  heat  feedback  to  the  propellant  surface  can  be  computed  from 
steady  state  burning  rate  data.  The  pressure  field  in  the  propellant 
bed  is  coupled  to  the  transient  heat  conduction  equation  in  the  burning 
solid  to  obtain  the  instantaneous  regression  rate.  Kooker  and  Nelson^ 
demonstrated  that  invariant  imbedding  can  be  used  to  solve  this  type 
combustion  model.  Gough^  provided  the  two-phase  flow  model  for  the  gas 
properties.  Coupling  the  two  requires  simultaneously  solving  for  the 
gas  pressure  and  the  regression  rate  of  the  solid. 


1.  E.  B.  Fisher  and  A . P.  Trippe , "Application  of  a Flame  Spread  Model 
to  Design  Problems  in  the  155mm  Propelling  Charge" , Proceedings  of 
the  12th  JANNAF  Combustion  Meeting , Newport , RI,  August  1975. 

2.  A.  W.  Horst , C.  W.  Nelson , and  I . W.  May,  "Flame  Spreading  in  Granu- 
lar Propellant  Beds:  A Diagnostic  Comparison  of  Theory  to  Experiment" , 
AIAA  Paper  77-856,  13th  AIAA  Propulsion  Specialists  Meeting,  Orlando, 
FL,  July  1977. 

3.  D.  E.  Kooker  and  C.  W.  Nelson,  "Numerical  Solution  of  Three  Solid 
Propellant  Combustion  Models  During  a Gun  Pressure  Transient" , Pro- 
ceedings of  the  12th  JANNAF  Combustion  Meeting,  Newport,  RI,  August 
1975. 

4.  L.  H.  Caveny,  M.  Summerfield,  and  C.  W.  Nelson,  "Ignition  Transients 
and  Pressurization  in  Closed  Chambers" , Proceedings  of  the  11th  JANNAF 
Combustion  Meeting,  Pasadena,  CA,  September  1974. 

5.  P.  S.  Gough,  "A  Quasi-One- Dimensional  Two  Phase  Flow  Model  of  Interior 
Ballistics:  Towards  a Universal  Gun  Code",  Proceedings  of  the  13th 
JANNAF  Combustion  Meeting,  Monterey,  CA,  September  1976. 
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II.  THEORY 


The  governing  equations  and  numerical  approach  for  the  hydrodynamics 
of  the  flow  have  been  described  earlier^>6#  of  interest  here  is  the 
coupling  of  the  transient  combustion  with  the  flow. 


In  the  flow  field  is  distributed  a set  cf  marker  particles,  fewer 
in  number  than  grid  points,  convected  by  the  propellant  grain  velocity. 
These  marker  particles  are  used  to  compute  the  transient  heating  to  ig- 
nition and  mass  generation  by  combustion  within  the  bed.  A two-way 
interpolation  operates  the  continual  shuttling  between  markers  and  grid 
points.  Conditions  at  the  flow  mesh  are  interpolated  to  compute  a 
heat  flux  to  the  marker  particle.  Once  the  surface  temperature  of  the 
particle  has  been  calculated,  the  temperatures  in  the  flow  field  are 
interpolated  for  the  points  between  bounding  marker  particles.  Con- 
sider a flow  mesh  with  points  at  x^,  x^>  x^,  and  marker  particles 

at  x and  x,  as  shown 
a b 

o o o o o o 

X1  xa  X2  X3  Xb  X4 

with  heat  flux  Q^,  Q0 , and  computed  from  flow  conditions.  Then  by 
interpolation 


■ % + «2  - 


(1) 


The  temperature  field  in  the  particle  is  then  computed  by  solving  the 

transient  heat  conduction  equation  to  obtain  surface  temperatures 

and  T,  . T0  is  then  obtained  again  by  interpolation 
b Z 


x - x 

T„  = T + (T,  - T )(— -) 

2 a b a x^  - x^ 


(2) 


Once  ignition  of  a bounding  particle  (x^  or  i^)  occurs,  the  procedure 
should  recognize  the  shift  from  a simple  convective  heat  transfer.  The 
flux  is  computed  as  though  the  solid  were  still  being  convectively 
heated  but  the  heat  is  not  taken  from  the  gas  phase  energy  equation. 
This  process  continues  until  the  entire  solid  phase  is  ignited.  In- 
accuracies admitted  by  this  procedure  are  accepted  to  avoid  the  compu- 
tational expense  of  tracking  the  temperature  field  in  the  entire  solid 
phase  region.  Typically  eight  marker  particles  were  used  in  a grid 
mesh  of  35  points. 


6.  P.  S.  Gough 3 "Numerical  Analysis  of  a Two-Phase  Flow  with  Explicit  In- 
ternal Boundaries" , Paul  Gough  Associates  Report  PGA-TR-76-2 j Septem- 
ber 1976  (See  also  NOS  Indian  Head  Report  IHCR  77-5 , April  1977). 
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Before  ignition,  the  solution  to  the  transient  heat  conduction 
equation  is  straightforward,  given  the  flux  to  the  heated  surface. 

After  ignition,  it  is  considerably  more  difficult  because  the  non-linear 
equation  has  a coupled  non-linear  boundary  condition.  For  the  Zeldovich 
transient  model  (4),  the  surface  gradient  is  given  by 


9T 

9Z 


= -[t  — - £n(— ) - t ] 


(3) 


The  iterative  procedure  assumes  a trial  value  for  the  regression  rate 
to  solve  for  the  surface  temperature  Ts  by  an  invariant  imbedding  inte- 
gration (also  iterative)  of  the  heat  conduction  equation.  The  com- 
puted value  of  T is  then  used  in  a global  burning  rate  expression 
s 

r = A exp [-E/RT  ] (4) 

s 

Successive  values  of  r are  computed  until  two  agree  within  a specified 
difference. 


Values  for  rs  in  Eq.  (3)  are  obtained  by  applying  the  quasi-steady 
power  law  (apn)  to  the  linearly  interpolated  pressure  at  x . 

9. 

This  assumption  of  a sudden  transition  at  the  ignition  temperature 
to  self-sustaining  combustion  contradicts  observed  propellant  ignition 
which  shows  a developing  flame  requiring  some  continued  external  stimu- 
lus  after  first  flame  appearance7*8.  Sustained  burning  was  observed 
to  be  preceded  by  a brief  flux-dependent  transient  flame  followed  by 
radiation  assisted  burning. 

There  is  little  direct  evidence  to  validate  a quantitative  transi- 
tion model  with  convective  heating  fluxes.  Qualitative  arguments, 
which  can  offer  mechanisms  for  the  transition  flux,  yield  no  numerical 
values  needed  for  a calculation.  NOVA's  combination  of  radiation  and 
forced  convection  necessitates  some  arbitrary  decisions  on  heat  feed- 
back during  the  transition. 


7.  T.  J.  Ohlemiller>  L.  H.  Caveny>  L . DeLuoa 3 and  M.  Summerfieldj 

rr Dynamic  Effects  on  Ignitability  Limits  of  Solid  Propellants  Sub- 
jected to  Radiative  Heating ",  Fourteenth  Symposium  (International) 
on  Combustion 3 The  Combustion  Institute 3 pp.  1297-1307 3 (1973). 

8.  L . DeLucaj  T.  J.  Ohlemiller > L . H . Caveny3  and  M.  Summer  fields 
"Radiative  Ignition  of  Double  Base  Fropellants:  II  Pre-Ignition 
Events  and  Source  Effects ",  AIAAJ , 14_  (8),  p.  1111-1117  (1976). 
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Andersen  suggests  one  criterion  for  the  time  to  self-sustaining 
combustion  as 


x = a/r 


2 

s 


(5) 


Using  this  criterion  requires  a value  of  rs  in  a transient  pressure  field. 
Estimates  of  the  local  pressurization  rate,  d?/dt,  can  be  used  to  extra- 
polate the  pressure.  Eq.  (5)  can  then  be  expressed  as 


T = 


a/ {a  [P 


ign 


2 


(6) 


After  a trial  and  error  solution  yields  the  transition  time,  a mostly 

arbitrary  decrease  in  convective  flux  augments  the  flame  heat  feedback. 

One  approximation  is  a linear  decrease  in  time  from  the  flux  at  ignition 

(Q.  ) such  that 

ign 


Q = Q.  (t  < t . + t) 

c xign  ign 


Total  heat  flux  to  the  surface  is  then 


(7) 


Q - A*  + Qc  • (8) 

After  the  transition  time  x has  expired,  heat  is  supplied  only  by  the 
i Zeldovich  flame. 

Kuo^  postulated  that  the  gas-to-particle  heat  transfer  coefficient 
decreased  linearly  with  surface  temperature, 


h = h 

o 


T. 

LT. 

ign 


(9) 


He  reported  a typical  transition  time  of  abouc  lys , which  is  far  smaller 
than  Andersen* s criterion  and  than  DeLuca’s®  observation  of  milliseconds. 


There  seems  little  immediate  prospect  of  obtaining  a definite  es- 
timate of  the  actual  heat  transfer  from  the  developing  flame.  KookerU 
promises  an  eventual  calculation  once  the  kinetics  can  be  identified. 


9.  W.  H.  Andersen "Model  of  Transient  Ignition  to  Self-Sustaining 
Combustion" Comb . Sci . and  Tech.  5_3  p.  ?53  1972 . 

10.  K.  K.  Kuo3  R.  Vichnevetsky  and  M.  Summerfield3  " Theory  of  Flame 
Front  Propagation  in  Porous  Propellant  Charges  under  Confinement"  3 
Princeton  Univ.  AMSE  Report  10003  August  1971  (AD  #762063). 

11.  D.  E.  Kooker3  "numerical  Predictions  of  Ignition  and  Flamespreading 
in  Ozone/Oxygen  Mixture" Proceedings  of  Ihe  14th  JANNAF  Combustion 
Meeting 3 Colorado  Spring s3  C03  August  1977. 
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As  will  be  seen  in  the  results,  the  Zeldovich  heat  feedback  by  itself 
is  not  enough  to  sustain  combustion  at  low  pressures. 

In  a typical  heating  to  ignition,  the  heating  rate  rises  £rom  zero 
to  640  J/cnr/sec  in  about  0.5ms.  The  average  rate  of  430  J/cm  /sec 
is  too  high  to  compare  directly  with  DeLuca's  experimental  results.  At 
DeLuca's  highest  rate  (80  cal/cm^/sec)  the  delay  time  between  faint  and 
strong  IR  signal  (Fig.  3a,  Ref.  8)  was  0.5ms  at  21atm  to  2ms  at  5atm. 

At  the  pressures  calculated  by  NOVA  where  ignition  occurs  at  about  30atm, 
the  transition  times  should  be  less  than  1 ms  to  a self-sustaining 
combustion. 

4 

Earlier  studies  by  Caveny  et  al  noted  the  transition  problem. 
Heating  was  imposed  at  840  J/cm^/sec  with  a linear  decay  to  zero  in 
0.5  or  1.0ms  after  ignition  which  occurred  at  0.5MPa.  To  insure  a 
continued  flame,  heat  from  a forced  convection  to  the  surface  was 
added  whenever  Zeldovich  stability  parameters  indicated  unstable 
burning.  Relative  burning  rate  peaked  at  a pressure  of  3.3MPa. 


Given  the  unavailability  of  a full  solution  and  the  arbitrary  pre- 
vious approaches,  it  is  not  unreasonable  to  adopt  the  Andersen  criterion 
with  a linear  ramp  as  given  by  Eq.  (7). 


A linear  interpolation  of  relative  burning  rate  (r/r  ) converts  the 
marker  particle  transient  rates  to  the  grid  mesh  rates 


r2  ■ + • 
2sa  sb  s a d a 


(10) 


Again  a judgment  is  made  that  relative  burning  rate  is  a more  useful  varia- 
ble to  interpolate  than  absolute  rate.  When  the  interpolation  region  in- 
cludes a steep  pressure  gradient,  neither  absolute  nor  relative  rates 
should  be  expected  to  vary  linearly  with  distance.  Although  it  is  not 
yet  clear  which  is  closer  to  linearity,  intuition  says  that  interpolating 
the  relative  rate  will  cause  a smaller  error.  A sensitivity  test  must 
precede  any  definitive  judgment. 


III.  UNCOUPLED  PROBLEM 


To  lower  computer  costs  while  investigating  the  transient  burning 
response  to  the  expected  gun  conditions,  the  NOVA  transient  burning  sub- 
routines were  extracted  and  run  with  imposed  typical  pressure  and  heat 
flux  histories.  Parametric  studies  then  established  the  sensitivity  of 
the  predicted  transient  burning  rate  response  to  variations  in  the  gun 
conditions  and  input  data.  For  the  155  with  nominal  input  data,  the 
heat  flux  computed  in  Ref.  12  and  pressure  at  the  marker  particle 
nearest  the  breech  could  be  approximated  as 


12.  C.  W.  Nelson > "Comparison  of  Predictions  of  Three  Two-Phase  Flow 
Codes"j  13th  JANNAF  Combustion  Mtg Monterey ^ CA^  September  1976 . 
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dP 

dt 


7000 


MPa 


(12) 


sec 

After  ignition,  heat  flux  is  computed  by  Eq.  (7)  and  Eq.  (8). 

Surprise!  The  propellant  flame  cannot  sustain  itself.  Heat  feed- 
back from  the  Zeldovich  formulation  is  not  enough  to  maintain  the  com- 
bustion unless  supplemented  by  an  assumed  source  outside  the  flame. 

The  delay  time  of  0.5ms  computed  from  Andersen’s  criterion  must  be 
arbitrarily  longer.  A series  of  trials  showed  that  about  2ms  was 
enough  to  allow  a transition  to  stable  behavior  in  this  problem.  Succes 
ful  transition  is  the  correct  criterion  because  the  gun  charge  does  ig- 
nite and  burn.  Such  pathological  examples  as  stuttering  or  delayed 
ignition  were  not  intended  to  be  handled  with  this  code  and  will  be 
left  for  another  study. 

Sensitivity  studies  were  performed  with  the  uncoupled  version  to 
find  the  combustion  guidelines  for  the  full  gun  cases.  With  the  155 
propellant  properties  of  Table  II  some  of  the  results  are  shown  in 
Table  I. 


Table  I. 

Sensitivity  Studies, 

Uncoupled  Problem 

Case 

Transition 
Time  (ms) 

OpCR)-1 

max  ( — ) 
r 

s 

Remarks 

i 

2.5 

0.004 

4.0 

2 

6.0 

0.003 

3.1 

3 

1.0 

0.003 

2.9 

4 

0.5 

0.003 

1.5 

Flameout 

5 

2.0 

0.002 

2.2 

6 

0.5 

0.002 

1.3 

Flameout 

7a 

6.0 

0.004 

4.0 

Flameout 

8 

6.0 

0.003 

2.9 

E 17  kcal/mol 

9 

2.0 

0.0025 

2.4 

s 

Note  a - Regression  rate  convergence  0.1  percent. 

It  is  discomforting  to  discover  a go-no  go  sensitivity  to  such  an 
arbitrary  parameter  as  the  details  of  an  indeterminable  heat  flux.  Be- 
cause the  combustion  details  are  secondary  to  the  hydrodynamics  of  the 
flow  in  a burning  propellant  bed,  the  transition  parameter  must  be  set 
to  avoid  extinguishment.  The  price  to  be  paid  is  a higher  peak  rela- 
tive burning  rate  and  a greater  mass  generation  in  the  transient  period. 


IV.  GUN  CALCULATIONS 


With  the  transient  burning  coupled  to  the  hydrodynamics,  NOVA  was 
used  to  calculate  the  performance  of  two  guns  for  which  tests  data  are 
available  for  comparison.  Earlier  reports  have  described  attempts  to 
simulate  the  gun  performance  by  quasi-steady  burning  rates  modified  (in 
some  cases)  to  obtain  agreement  with  test  data.  Input  data  for  the  two 
guns  are  shown  in  Tables  II  and  III. 


Table  II.  Input  Data 


Initial  Temperature 
Propellant  Density 
Propellant  Charge 
Grain  Shape 
Outer  Diameter 
Perforation  Diameter 
Grain  Length 

Chemical  Energy  of  Propellant 
Gas  Molecular  Weight 
Specific  Heat  Ratio 
Co-volume 

Igniter  Chemical  Energy 
Igniter  Molecular  Weight 
Igniter  Specific  Heat  Ratio 
Tube  Diameter  at  Breech 
Tube  Radius  at  Bore 
Projectile  Mass 
Settling  Porosity  of  Bed 
Bed  Rear  Boundary 
Bed  Forward  Boundary 
Primer  Input : 

Length 
Flow  Rate 


155mm  Howitzer 

294  K 

1.58  g/cm 
9.9  kg 

7 perf  cylinder 
10.67  mm 
0.86  mm 
24.3  ram 
4420  J/g 
23.46 
1.24 

.944  cm  /g 
. 3980  J/g 

23.0 
1.245 

185  mm 

155  mm 

43.1  kg 
0.53 
6.4  mm 

553  mm 


5/54  Gun 

306  K 

1.55  g/cm 
9.4  kg 

7 perf  cylinder 

10.0  mm 
1. 0 mm 

23.4  mm 
3900  J/g 

22.89 
1.246 
.933  cm/ g 
1568  J/g 

36.1 
1.25 

142  mm 

127  ram 

24.1  kg 
0.50 
0. 0 mm 

807  mm 


51mm  51mm 

.27  kg/mm/sec  .05  kg/mm/sec 
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Table  III.  Propellant  Combustion  Parameters 


Activation  Energy  (KJ/mol) 
Pre-exponential  (cm/ sec) 

Reference  Surf  Temp  (K) 

Thermal  Conductivity  (J/cm  /secK) 
Thermal  Diffusivity  (cm^/sec) 
Ignition  Temperature  (K) 
Temperature  Sensitivity  ( /K) 
Steady  State  Burning  Rate 
Coefficient  cm/s(MPa)n 
Exponent 


155mm  Howitzer  5/54  Gun 


62.8 

62.8 

2.43  x 

io5 

7.07  x 

10‘ 

623 

—A 

650 

8.73  x 

10-4 

8.73  x 

10 

8.68  x 

10  4 

8.68  x 

10 

450 

428 

0.0054 

0.0045 

0.387 

0.0598 

0.7 

1.2294 

0.597 

1 

0.4194 

1 

0.178 

0.7965  C 


a - P < 17.2  MPa 

b - 17.2  < P < 37.9  MPa 
c - 37.9  < P < 690  MPa 

V.  155mm  HOWITZER 

The  standard  155mm  howitzer  charge  was  made  into  a full  bore  bag, 
an  experimental  charge  to  approach  one  dimensizmal  geometry.  No  free 
space  exists  between  bag  and  chamber  wall.  For  the  given  propellant 
weight,  the  charge  was  shorter  (55cm)  than  the  standard  bag  charge 
(76cm).  The  only  other  change  from  the  standard  charge  was  the  removal 
of  the  center  core  igniter  to  permit  only  base  ignition. 

Earlier  NOVA  calculations  of  the  155  simulated  the  gun  pressure 
wave  behavior  when  a constant  was  added  to  the  quasi-steady  regression 
rate  (as  it  was  in  a similar  5/54  gm  calculation^).  A constant  0.5  cm/ 
sec  gave  an  acceptable  match  with  a typical  experimental  record.  At 
low  pressure  (less  than  5 MPa)  the  regression  rate  then  exceeded  four 
times  the  quasi-steady  rate.  At  pressures  above  50  MPa,  the  multiple 
was  less  than  1.5  which  is  smaller  but  still  not  insignificant.  The 
effects  of  the  increase  are  felt  in  the  first  cycle  of  the  wave; 
damping  thereafter  is  qualitatively  the  same  as  in  the  quasi-steady 
rate  case. 
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Later  examination  of  the  additive  constant  simulations  showed  that 
adjustment  of  the  bore  resistance  and  projectile  seating  distance  were 
also  sufficient  for  an  adequate  simulation.  Projectile  seating  dis- 
tance is  a measurable  input  and  should  not  be  an  arbitrary  variable. 

It  was  unfortunate  that  it  was  incorrectly  stated  in  some  simulations. 
What  resulted  was  the  finding  that  it  could  act  as  an  adjustable  con- 
stant. An  additive  burning  rate  constant  had  a smaller  effect  on  the 
pressure  wave  amplitudes.  Figure  1 shows  the  simulations  for  the  re- 
vised input  data  [curve  labeled  (QS ) ] . Shaded  portion  shows  test  re- 
sults. Amplitudes  for  the  earlier  version  with  an  additive  constant  of 
1.27  cm/sec  were  only  marginally  (10%)  higher.  An  unpleasant  conclusion 
drawn  from  this  comparison  is  that  either  projectile  motion  or  regression 
rate  may  be  arbitrarily  adjusted  to  force  simulations  to  match  test 
data. 


Bore  resistance  is  derived  from  Picatinny  Arsenal  data  for  155mm 
howitzer  with  charges  zone  3 and  4.  The  data  scatter,  especially  at  low 
projectile  velocities,  causes  a relatively  wide  interval  for  a resistance 
estimate.  The  actual  values  used  in  a linear  interpolation  are  as 
follows : 

Distance  Traveled  (cm)  0 1.5  9.1  514 

Resistance  Pressure  (MPa)  6.9  41.4  17.2  10.3 

The  principal  effect  of  the  transient  regression  was  to  increase  the 
amplitude  of  pressure  wave  beyond  the  already  too  high  predictions  with 
the  quasi-steady  regression.  Figure  1 compares  predictions  with  experi- 
ment (two  tests).  Flame  speed  of  about  0.4mm/ys  was  unaffected 
(Figure  2). 

One  variation  changed  the  quasi-steady  burning  rate  dependence  by 
increasing  the  pressure  exponent  to  0.85  while  retaining  the  same 
Arrhenius  activation  energy  and  reference  surface  temperature.  Earlier 
calculations  by  Horst^  in  a 5 inch  gun  predicted  a lower  amplitude 
with  higher  exponent  because  the  effect  is  to  depress  the  low  pressure 
regression  rates  and  to  increase  the  burning  rate  at  pressures  above 
that  where  the  rate  is  unchanged,  usually  called  the  pivot  pressure. 

(With  this  propellant,  the  burning  rate  at  7MPa  drops  from  1.5  cm/ sec  to 
.72  cm/sec.)  Pressure  wave  amplitudes  decreased  as  seen  in  Figure  1. 


13.  J.  DeLorenzo  and  A.  C . Vallado 3 " Compilation  of  Traces  and 
Tabulation  of  Round- to- Round  Data  for  the  155mm  XM211  Program 
(Phase  I)  Conducted  at  Picatinny  Arsenal  12-13  April  19763 
unpublished 3 Picatinny  Arsenal 3 May  1976 . 

14.  A.  W.  Horst 3 n Influence  of  Propellant  Burning  Rate  Representation 
on  Gun  Environment  Flamespread  and  Pressure  Wave  Predictions ", 
IHMR  76-255 3 March  1976. 
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Figure  1.  Pressure  Wave  Predictions  for  155mm  Howitzer 
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DISTANCE  FROM  BREECH  (cm) 

Figure  2.  Flame  Spreading  Predictions 


VI.  5/54  GUN 


The  standard  Navy  5/54  propelling  charge  and  projectile  were  modi- 
fied by  (1)  replacing  the  packing  elements  between  charge  and  projectile 
by  an  incompressible  plastic  disc  to  eliminate  effects  of  compressible 
packing  and  (2)  replacing  the  metallic  rotating  band  with  a plastic 
band  to  eliminate  almost  all  rotating  band  resistance  to  motion.  Two 
igniter  configurations  were  fired:  the  standard  bayonet  igniter,  and 
a black  powder  base  bad  igniter. 

(A  note  on  nomenclature.  "Igniter"  here  refers  to  the  black  powder 
or  other  material  which  provides  the  hot  gases  to  ignite  the  propellant 
grains.  There  is  a difference  in  Army  and  Navy  use  of  "primer"  and 
"igniter".) 

Earlier  calculations  of  the  base  pad  igniter  with  an  additive  burning 
constant  of  1.27  cm/sec,  successfully  simulated  the  test  behavior  (2). 

It  should  be  noted  here  that  the  transient  burning  model  used  in 
the  5/54  calculations  differed  from  that  in  the  155  calculations  in  that 
there  is  no  transition  to  a self-sustaining  flame.  Convective  heat 
transfer  from  the  chamber  gases  to  the  propellant  surface  is  assumed  to 
continue  through  the  propellant  flame.  This  heat  is  added  to  the  flame 
heat  feedback  as  computed  by  the  Zeldovich  model.  Regression  rates  are 
thus  higher  than  the  quasi-steady  rates  even  when  any  pressure  transients 
end.  This  treatment  gives  a limiting  analysis  to  show  the  worst  case 
of  combined  transient  and  erosive  burning. 

Although  it  avoids  an  arbitrary  decision  on  flame  development,  the 
results  show  it  to  be  "too  worst"  to  the  extent  that  credible  flow  calcu- 
lations are  not  obtained.  The  calculated  pressure  waves  exceed  the  code’s 
built-in  limitations  and  thus  the  calculations  stop  on  steep  pressure 
transients  during  bed  stagnation  at  the  projectile  base.  Figure  3 shows 
the  wave  development  for  both  igniters. 

Explanation  can  be  found  in  the  relative  regression  rates  of  Figure 
4.  For  the  155  howitzer  the  relative  rate  peaks  around  two  and  decays 
to  the  quasi-steady  rate  (relative  rate  1.0).  For  the  base  ignited  5/54 
gun,  the  relative  rate  rises  to  around  five  and  then  oscillates  as  it 
decays  toward  the  quasi-steady  rate.  This  oscillation  correlates  directly 
with  the  parameter  |u-Up|  which  says  that  the  convection  from  the  chamber 
gas  strongly  contributes  to  the  heat  feedback  to  the  propellant  surface. 
Such  an  overstated  regression  rate  creates  a stronger  pressure  front 
which  compacts  the  grains  harder  against  the  projectile  base. 

The  effect  on  the  flame  spreading  can  be  seen  in  Figure  2.  The 
quasi-steady  regression  and  the  cubic  profile  method  for  heating  to  ig- 
nition combine  for  a flame  spreading  rate  of  around  0.25mm/ys,  not  far 
from  the  sound  speed  in  the  undisturbed  gas.  The  fully  transient  solu- 
tion produces  a slow  initial  spreading  (.04mm/ys)  at  the  breech  end 
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Figure  3.  Pressure  Differences  for  5/54  Gun 
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accelerating  to  0.4  mm/ys  at  the  forward  end.  The  cause  for  the  differ- 
ence is  not  clear.  A major  difference  between  the  two  guns  is  that  the 
5/54  has  a lower  igniter  input  rate.  How  this  fact  couples  with  the 
heat  conduction  solution  remains  to  be  determined.  A summary  of  effects 
for  both  guns  can  be  seen  in  Table  IV. 


Table  IV.  Summary  of  Results 


Burning 

Flame 

AP* 

p 

Case 

Rate 

Igniter 

Speed  (mm/ys) 

(MPa) 

max (MPa) 

155mm  Howitzer 

i 

Steady 

Base 

0.4 

- 73/39 

293 

IT 

Transient 

Base 

0.4 

-90/76 

NC 

2 

Add  0 . 5 cm/  s 

Base 

0.4-0. 5 

-83/56 

308 

Experimental 

Base 

Unknown 

-40/35 

235 

3 

n = 0.85 

Base 

0.3 

-50/NC 

NC 

5/54  Gun 

11 

Steady 

Base 

0.1-0.25 

-24/40 

235 

11T 

Transient 

Base 

0.04-.4 

NC 

NC 

12 

Steady 

Bayonet (22")  0.4 

-8/17 

241 

12T 

Transient 

Bayonet 

(22")  0.4 

NC 

NC 

* Breech  Pressure  - Forward  Pressure.  First  number  is  maximum 
first  negative  difference;  second  number  is  maximum  second  positive 
difference. 

The  bayonet  primer  configuration  exhibits  no  such  differences.  As 
seen  in  Figure  5,  the  rates  are  all  about  the  same  and  equal  to  the  rate 
measured  by  McClure  and  East-^  with  M26  propellant.  (Such  agreement 
may  be  fortuitous  given  the  differences  in  grain  sizes  and  composition). 

VII.  THE  CONVERGENCE  BARRIER 

The  non-linear  equation  and  boundary  condition  require  an  iteration 
to  obtain  a consistent  set  for  burning  rate,  surface  temperature,  and 
heat  flux.  Kooker  and  Nelson  used  a convergence  criterion  of  about  10  K 


15.  D . /?.  McClure  and  J . L . East , " Experimental  Techniques  for  Investi- 
gating the  Start-Up  Ignition/ Combustion  Transients  in  Full  Scale 
Charge  Assemblies" 3 11th  JANNAF  Combustion  Meeting 3 Pasadena 3 CA3 
September  1974 . 
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for  surface  temperature  and  10  cm/sec  for  burning  rate.  Later  study 
showed  that  relaxing  these  by  two  orders  of  magnitude  caused  a 10%  change 
in  peak  relative  burning  rate  but,  more  usefully,  only  a 1%  error  in  gas 
generation  during  the  transient.  The  earlier  Zeldovich  study  by  Caveny 
et  al^  used  an  explicit  scheme  without  iteration  by  using  the  last 
computed  value  for  the  regression  rate. 

-4 

At  first,  the  convergence  criteria  were  set  at  10  percent  of 
regression  rate  and  surface  temperature.  These  translate  to  about 
10“^K  and  10-~*  cm/sec  during  the  transients.  It  was  soon  discovered 
that  convergence  was  not  obtained  in  ten  iterations  when  the  regression 
rate  exceeded  about  3 cm/ sec. 

Because  it  was  too  costly  to  decrease  the  problem  time  step  to 
accommodate  the  regression  transient,  an  operator  splitting  technique 
was  adopted^  as  had  been  previously  used  by  Dwyerl?.  Each  hydro- 
dynamic  time  step  is  split  into  the  necessary  sub-steps  to  retain  con- 
vergence of  the  regression  rate  solution  with  a MacCormack  scheme  Lf 
for  the  hydrodynamics  and  a Gear  Scheme  Lc  for  the  chemistry, 

Dwyerl7  applied  the  operators  in  a sequence 


n+2  T T T t n 
u . = L _L  L L,_u. 

l f c c f 1 


(13) 


to  advance  the  dependent  variable  u two  time  steps.  This  present  study 
used  the  simple  scheme 


n+1  T T n 
u.  = L _L  u.  3 
l f c 1 


(14) 


which  is  a crude  approximation  of  a scheme  originally  intended  for  multi- 
dimensional problems.  Some  improvement  can  be  made  by  interpolating  the 
pressure  field  across  the  interval  for  each  time  sub-step.  Other  users 
have  assigned  a single  representative  value  of  the  hydrodynamic  variables 
for  all  the  sub-steps.  The  Navy  gun  calculations  used  the  interpolated 
values;  the  Army  calculations  used  only  the  forward  value.  Sensitivity 
to  the  use  of  either  values  has  yet  to  be  determined. 


16.  N.  N.  Yanenko,  The  Method  of  Fractional  Steps , translated  by 
M . Holt,  Spring er-Ver lag,  Berlin,  1971 . 

17.  H.  A.  Dwyer  and  B.  R.  Sanders,  "Modeling  of  Unsteady  Combustion 
Phenomena ",  AIAA  Paper  77-136,  15th  Aerospace  Sciences  Meeting, 
January  1977. 
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Even  the  time-splitting  technique  was  inadequate  for  the  task  when 
the  regression  rate  exceeded  about  10  cm/sec.  The  same  trouble  was 
noted  by  Kooker  and  Nelson  in  that  the  time  step  had  to  be  cut  to 
about  10~8  sec  when  the  rate  was  about  6 cm/sec.  It  is  clear  that 
some  drastic  changes  are  needed  if  accurate  and  affordable  calculations 
of  the  transient  rate  are  to  be  made  at  high  regression  rates.  For 
the  present  problem,  the  convergence  failure  occurs  as  the  compression 
wave  strikes  the  projectile  base  and  the  pressure  nears  150  MPa  where 
the  quasi-steady  regression  rate  is  about  12  cm/sec  and  the  transient 
rate  about  13  cm/sec.  There  a time  step  of  10”  ^ sec  is  too  large. 

There  can  be  no  temporary  acceptance  of  the  small  time  step  at  the 
local  high  pressure  region  because  the  whole  bed  will  soon  exceed  200 
MPa.  A compromise  can  recognize  that  for  quasi-steady  regression  rates 
above  10  cm/ sec,  the  transients  may  be  ignored.  Previous  studies^  pre- 
dicted that  the  transient  rate  approaches  the  quasi-steady  rate  for 
pressures  above  about  50  MPa. 

Thus,  those  transient  calculations  which  go  to  completion  have 
had  the  transient  combustion  terminated  whenever  the  quasi-steady 
regression  rate  exceeded  10  cm/sec.  If  a transient  still  exists  at 
that  pressure  (or  will  be  later  created  by  the  passage  of  a steep  com- 
pression wave  through  the  bed)  an  error  is  thereby  introduced.  For 
now  it  is  assumed  that  the  error  can  be  neglected. 

A new  iteration  algorithm  is  needed  before  fully  transient  calcu- 
lations can  be  completed. 


VIII.  CONCLUSIONS 

1.  Transient  burning  can  affect  interior  ballistics  calculations. 

2.  Zeldovich  model  response  agrees  with  uncoupled  response  in  similar 
pressure  history. 

3.  For  practical  IB  calculations,  cubic  profile  may  be  unsatisfactory  for 
flame  spreading  with  weak  igniters. 

4.  Further  evaluation  awaits  improvements  in  convergence  and  transition 
to  self-sustaining  combustion. 


18.  C.  W.  Nelson^  "Response  of  Three  Types  of  Transient  Combustion 
Models  to  Gun  Pressurization" 3 BRL  Memorandum  Report  No.  27 523 
May  1977 3 AD  #A041079. 
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